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We study the entire coupled evolution ol the inflaton <f)(t) and the scale lactor a(t) for general 
initial conditions <j>(to) and d<j>(to)/dt at a given initial time to. The generic early universe evolution 
has three stages: decelerated fast-roll followed by inflationary fast-roll and then inflationary slow- 
roll (an attractor always reached for generic initial conditions). This evolution is valid for all 
regular inflaton potentials v(4>). In addition, we find a special (extreme) slow-roll solution starting 
at t — — oo in which the fast-roll stages are absent. At some time t = i«, the evolution backwards 
in time from to reaches generically a mathematical singularity where a(t) vanishes and the Hubble 
parameter becomes singular. We determine the general behaviour near the singularity. The classical 
homogeneous inflaton description turns to be valid for t — t, > 10 tpi anc k well before the beginning 
of inflation, quantum loop effects are negligible there. The singularity is never reached in the 
validity region of the classical treatment and therefore it is not a real physical phenomenon 
here. Fast-roll and slow-roll regimes are analyzed in detail including the equation of state evolution, 
both analytically and numerically. The characteristic time scale of the fast-roll era turns to be 
t\ = (1/m) -\/V(0)/[3 M 4 ] ~ 10 4 t pianck where V is the double- well inflaton potential, m is the 
inflaton mass and M the energy scale of inflation. The whole evolution of the fluctuations along 
the decelerated and inflationary fast-roll and slow-roll eras is computed. The Bunch-Davies initial 
conditions (BDic) are generalized for the present case in which the potential felt by the fluctuations 
can never be neglected. The fluctuations feel a singular attractive potential near the t = t, 
singularity (as in the case of a particle in a central singular potential) with exactly the critical 
strength (—1/4) allowing the fall to the centre. Precisely, the fluctuations exhibit logarithmic 
behaviour describing the fall to t = t*. The power spectrum gets dynamically modified by the 
effect of the fast-roll eras and the choice of BDic at a finite time through the transfer function 
D(k) of initial conditions. The power spectrum vanishes at k = 0. D(k) presents a first peak 
for k ~ 2/770 (?7o being the conformal initial time), then oscillates with decreasing amplitude and 
vanishes asymptotically for k — > 00. The transfer function D(k) affects the low CMB multipoles Ct- 
the change ACi/Ce for 1 < £ < 5 is computed as a function of the starting instant of the fluctuations 
to. CMB quadrupole observations indicate large suppressions which are well reproduced for the 
range to — t, > 0.05/m ~ 10100 t P i anck . 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



Since the Universe expands exponentially fast during inflation, gradients are exponentially erased and can be 
neglected. At the same time, the exponential stretching of spatial lengths classicalizes the physics and allows a 
classical treatment. One can therefore consider a homogeneous and classical inflaton field which thus determines 
self-consistently a homogenous and isotropic Friedman-Robertson Walker metric sourced by this inflaton. 

This treatment is valid for early times well after the Planck time t = 10 -44 sec, at which the quantum fluctuations 
are expected to be large and thus a full quantum gravity treatment is required. 

In this paper we study the entire coupled evolution of the inflaton field ip(t) and the scale factor a(t) of the metric 
for generic initial conditions, fixed by the values of 4>{to) and d(f>(to)/dt at a given initial time to- 

We show that the generic early universe evolution has three stages: a decelerated fast-roll stage followed by an 
inflationary fast-roll stage and then by a slow-roll inflationary regime which is an attractor always reached for generic 
initial conditions. This evolution is valid for all regular inflaton potentials. In addition, we find a particular (extreme) 
slow-roll solution starting from t = — oo in which the fast-roll stages are absent. 

The evolution backwards in time from to reachs generically a mathematical singularity at some time t — t* where 
the scale factor a(t) vanishes, and the Hubble parameter becomes singular. 

We find the general behaviour of the inflaton and the scale factor near the singularity as given by eqs. (|2.14j) - (|2.17[) 
and determine the validity of the classical approximation, namely (H/Mpi) 2 <C 1. It must be stressed that such 
mathematical singularity is attained extrapolating the classical treatment where it is no more valid. The singularity 
is never reached in the validity region of the classical treatment and therefore such mathematical singularity is not 
a real physical phenomenon here. 

Quantum loops effects turns to be less than 1% for t — t* > 10~ 42 sec and therefore the classical treatment of the 
inflaton and the space-time can be trusted well before the begining of inflation. 

The fast-roll (both decelerated and inflationary) and slow-roll regimes are analyzed in detail, with both the exact 
numerical evolution and an analytic approximation, and the whole equation of state evolution in the three regimes. 
We consider here the double well (broken symmetric) fourth order inflaton potential since it gives the best description 
of the CMB+LSS data 2, 7] within the Ginsburg-Landau effective theory approach we follow. 

The characteristic time scale of the fast-roll era turns to be t\ — (1/m) y/V(0)/[3 M A ] ~ 10 4 tpianck where V(0) 
is the double well inflaton potential at zero inflaton field, m is the inflaton mass and M the energy scale of inflation. 
The time scale of the inflaton in the extreme slow roll solution goes as the inverse of ti, namely l/[m 2 t\]. 

We study the whole evolution of the curvature and tensor fluctuations along the three succesive regimes: decelerated 
fast-roll followed by inflationary fast-roll and then inflationary slow-roll, and compute the power spectrum by the end 
of inflation. The fluctuations feel a singular attractive potential near the t = singularity (as in the case of a 
particle in a central singular potential) with exactly the critical strength (—1/4) for which the fall to the centre 



3 



becomes possible. Precisely, the logarithmic behaviour of the fluctuations for t — > t* eq. (|4.3l) describes the fall to 
t = t* for the critical strength of the potential W-r felt by the fluctuations. 

We generalize the Bunch-Davies initial conditions (BDic) to the present case in which the potential felt by the 
fluctuations can never be neglected. 

In general, the mode functions for large k behave as free modes since the potential Wr becomes negligible in this 
limit except at the singularity t = t*. One can then impose Bunch-Davies conditions for large k which corresponds 
to assume an initial quantum vacuum Fock state, empty of curvature excitations 

\l 1 k 



and therefore 

dSn 
drj 



(k; rjo) k -i k S n (k; rjo) 



Here rj stands for the conformal time: dr\ = dt/a(t). Eq. (|l.ll) fulfils the Wronskian normalization (that ensures the 
canonical commutation relations) 

W[S n ,S* n ] = Sn^-^S^ = i. (1.2) 

In asymptotically flat (or conformally flat) regions of the space-time the potential felt by the fluctuations Wr(t]) 
vanishes and the fluctuations exhibit a plane wave behaviour for all k (not necesarily large). This is not the case in 
strong gravity fields or near curvature singularities as in the present cosmological space-time where W-r. (rj) can never 
be neglected at fixed k. However, we can choose Bunch-Davies initial conditions (BDic) at r\ = r] (or equivalently, 
t = to) by imposing 

dS 

— — (k;r]o) = -i k S-R(k;T]o) for all k . (1.3) 
drj 

That is, we consider the initial value problem for the mode functions giving the values of S-R(k;r]) and dS-R/drj at 
rj = rjo. This condition combined with the Wronskian condition eq. (jl.2l) implies that 



\Sn(k;r] )\ 



dS-R 

(«; vo, 



dr\ 



\ ■ (1-4) 



'I k 

which is equivalent to cq. (jl.ll) for large k. 

The power spectrum at the end of slow-roll inflation P-R(k) gets dynamically modified by the effect of the preceding 
fast-roll eras through the transfer function of initial conditions D{k): 

Pn(k) = P£ D (k)[l + D(k)] , (1.5) 

D(k) accounts for the effect of both the initial conditions and the fluctuations evolution during fast- roll (before 
slow-roll). D(k) depends on the time to at which BDic are imposed. 

The power spectrum P^ D (k) corresponds to start the evolution with pure slow-roll from to — > —oo and with BDic 
cq. ([1.3[) -cq. (ll.4[) imposed there at t — > — oo, that is 770 = —00. P-R D {k) is given by its customary pure slow-roll 
expression, 

~)BD / / \ _i„„ a n, \ 1 1\ i„„ k , 1 „ 1„„2 k , /o ( 1 



lo gj P^(fc) = log A s {ko) + (n s - 1) log — + i n run log z — +0 . (1.6) 

where N is the number of inflation cfolds since the pivot CMB scale ko exits the horizon. We take here N = 60. 

Actually, BDic can be imposed at rj = rjo = —00 if and only if the inflaton evolution also starts at rj = rjo = —00. 
This only happens for a particular inflaton solution: the extreme slow-roll solution that we explicitly present and 
analyze in sec. MI Al In the extreme slow-roll case the fast-roll eras are absent, BDic are imposed at to —> —00 (that 
is ?/o = ~ 00)7 then D(k) = and P-R(k) — P 7 f Z5 (/c). Only in this case the fluctuation power spectrum at the end of 
inflation is the usual power spectrum P^ D (k) eq. (|4.24j) . 
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When BDic are imposed at finite times to, the spectrum is not the usual P^ D {k) but it gets modified by a 
non-zero transfer function D(k) eq. (|4.2f \ . The power spectrum Pn(k) vanishes at k = and exhibits oscillations 
which vanish at large k [see figs. Eland [7] 

Generically, the power spectrum vanishes at k — and we thus have 

f + D(k) k =° 0{k n ° +1 ) . (1.7) 

as shown in sec. IV Al D(k) presents a first peak for k ~ 2/7/0 and then oscillates asymptotically with decreasing 
amplitude such that 

#) fc =°°o(^- (1.8) 



We solved numerically the fluctuations equation with the BDic eq. (|4.7[) covering both the fast- roll and slow-roll 
regimes, namely for different initial times to ranging from the singularity r = t* till the transition time T trans from 
fast-roll to slow-roll. That is to say, we solved the fluctuations evolution for BDic imposed at different times in 
the three eras and we compare the resulting power spectra among them. We computed the corresponding transfer 
function, D(k) for the BDic imposed at the different eras. We depict 1 + D(k) vs. k for the different values of the 
time to where BDic are imposed in figs. El 

When the BDic are imposed during the fast-roll stage well before it ends, D(k) changes much more significantly 
than along the extreme slow roll solution. This is due to two main effects: the potential felt by the fluctuations is 
attractive during fast-roll, and 770, (far from being almost proportional to l/a(r])), tends to the constant value 77* as 
t — > t + and a(rf) — > 0. The numerical transfer functions 1 + D(k) obtained from eqs. (|4.12[) and (14.211) are plotted in 

figs. El 

We have also computed D(k) analytically with BDic at finite times 770, and a simple form is obtained in the scale 
invariant case, which is the leading term in the slow-roll expansion: 

cos2x sin 2a; sin 2 ie , , _ 

D(k) = — x=kr,o. 1.9 

x z x i 



Different initial times to lead essentially to a rescaling of k in D(k) by a factor 770 since the conformal time 77 is 
almost proportional to 1/0(77) during slow-roll [see figs. Eland below eq. (|5.7|) ] . By virtue of the dynamical attractor 
character of slow-roll, the power spectrum when the BDic are imposed at a finite time to cannot really distinguish 
between the extreme slow-roll solution or any other solution which is attracted to slow-roll well before the time to- 

Using the transfer function D(k) we obtained, we computed the change on the CMB multipoles ACi/Ce for £ = 1,2 
and 3 as functions of the starting instant of the fluctuations to- We plot ACg/Ce for 1 < £ < 5 vs. to — t* in fig. 
|H We see that ACn/Ci is positive for small to — t* and decreases with to becoming then negative. The CMB 
quadrupole observations indicate a large suppression thus indicating that to — i* > 0.05/m ~ 10100 ip; a ncfe- 

The fact that choosing BDic leads to a primordial power and its respective CMB multipoles which correctly 
reproduce the observed spectrum justifies the use of BDic. 

Besides finding a CMB quadrupole suppression in agreement with observations we provide here predictions 

for the dipole and £ < 5-multipole suppressions. Forthcoming CMB observations can provide better data to confront 
our CMB multipole suppression predictions. It will be extremely interesting to measure the primordial dipole and 
compare with our predicted value. 



II. THE PRE- INFLATIONARY AND INFLATIONARY FAST-ROLL ERAS 



The current WMAP data are validating the single field slow-roll scenario [l]. Single field slow-roll models provide 
an appealing, simple and fairly generic description of inflation. This inflationary scenario can be implemented using 
a scalar field, the inflaton with a Lagrangian density (see for example ref. @) 



C = a 3 (t) 



2 



(Vy) 2 

2a 2 (t) 



V(<p) 



(2.1) 



5 



where V(ip) is the inflaton potential. Since the universe expands exponentially fast during inflation, gradient terms 
are exponentially suppressed and can be neglected. At the same time, the exponential stretching of spatial lengths 
classicalize the physics and permits a classical treatment. One can therefore consider an homogeneous and classical 
inflaton field tp(t) which obeys the evolution equation 

tp + 3H(t) ip + V'((p) =0 (2.2) 

in the isotropic and homogeneous FRW metric which is sourced by the inflaton 

ds 2 = dt 2 - a 2 (t) dx 2 (2.3) 

H(t) = a(t)/a(t) stands for the Hubble parameter. The energy density and the pressure for a spatially homogeneous 
inflaton are given by 

. 2 -2 

P=y+%) , P=\-V{V). (2.4) 



Thrcfore, the scale factor a(t) obeys the Friedmann equation, 

1 



H 2 (t) = — L- 



-^ 2 + V(<p) 



(2.5) 



In order to have a finite number of inflation efolds, the inflaton potential V((p) must vanish at its absolute minimum 

V'(cp mm ) = V(cp mm ) = (2.6) 

Otherwise, inflation continues forever. 

We formulate inflation as an effective field theory within the Ginsburg-Landau spirit 0, E3 • The theory of the 
second order phase transitions, the Ginsburg-Landau theory of superconductivity, the current-current Fermi theory 
of weak interactions, the sigma model of pions, nucleons (as skyrmions) and photons are all successful effective field 
theories. Our work shows how powerful is the effective theory of inflation to predict observable quantities that 
can be or will be soon contrasted with experiments. 

The effective theory of inflation should be the low energy limit of a microscopic fundamental theory not yet 
precisely known. The energy scale of inflation M should be at the Grand Unified Theory (GUT) energy scale in order 
to reproduce the amplitude of the CMB anisotropics [2J. Therefore, the microscopic theory of inflation is expected 
to be a GUT in a cosmological space-time. Such a theory of inflation would contain many fields of various spins. 
However, in order to have a homogeneous and isotropic universe the expectation value of the energy-momentum tensor 
of the fields must be homogeneous and isotropic. The inflaton field in the effective theory may be a coarse-grained 
average of fundamental scalar fields, or a composite (bound state) of fundamental fields of higher spin, just as in 
superconductivity. The inflaton does not need to be a fundamental field, for example it may emerge as a condensate 
of fermion-antifermion pairs < > in a GUT in the cosmological background. In order to describe the cosmological 
evolution is enough to consider the effective dynamics of such condensates. The relation between the effective field 
theory of inflation and the microscopic fundamental GUT is akin to the relation between the effective Ginzburg- 
Landau theory of superconductivity and the microscopic BCS theory, or like the relation of the 0(4) sigma model, 
an effective low energy theory of pions, photons and chiral condensates with quantum chromodynamics (QCD) [l6j . 

Vector fields have been considered to describe inflation in ref. 14]. The results for the inflaton should not be very 
different from the effective inflaton description since the energy-momentum tensor of the vector field is to be taken 
homogeneous and isotropic. Namely, we are always in the presence of a scalar condensate. 

Since the mass of the inflaton is given by M 2 /Mpi ~ 10 13 GeV massless fields alone cannot describe inflation 
which leads to the observed amplitude of the CMB anisotropies. 

The classical inflaton potential V(ip) gets modified by quantum loop corrections. We computed relevant quantum 
loop corrections to inflationary dynamics in ref. [2|, 1 1 51 ] . A thorough study of the effect of quantum fluctuations reveals 
that these l oop corrections are suppressed by powers of (H/Mpi) 2 ~ 10~ 9 where H is the Hubble parameter during 
inflation @, [H[. Therefore, quantum loop corrections are very small, a conclusion that validates the reliability of the 
classical approximation and of the effective field theory approach to inflationary dynamics. In particular, the (small) 
one-loop corrections to the potential in an inflationary universe are very different from the Coleman- Weinberg form 
id. 

We choose the inflaton field initially homogeneous which ensures it is always homogeneous. The fluctuations around 
are small and give small corrections to the homogeneity of the Universe. The rapid expansion of the Universe, in 
the inflationary regimes, takes care of the classical fluctuations, quickly flattening an eventually non-homogeneous 
condensate. 
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A. The complete inflaton evolution through the different eras 



It is convenient to use the dimensionless variables to analyze the inflaton evolution equations eas. (|2.2p ~ (I2.5|) . 



m 



t = mt 

The inflaton potential has then the universal form 

V(<p) = M 4 v 



M 



pi 



(2.7) 



Mpi 



(2.8) 



where M is the energy scale of inflation and v((f>) is a dimensionless function. Without loss of generality we can set 
v'(0) = Moreover, provided V"(0) ^ we can set without loss of generality |v"(0)| = 1/2. Namely, we have for 
small fields, 



;(0) ^^T-^+O^ 3 ) 



(2.9) 



where the minus sign in the quadratic term corresponds to new inflation and the plus sign to chaotic inflation. 

In these dimensionless variables, the energy density and the pressure for a spatially homogeneous inflaton are given 
from eq. (|2.4|) by 



M 4 2 1 dr 



■v{<f>) , 



P 



1 



v{ct>) , 



M 4 2 \dr y 

and the coupled inflaton evolution equation ()2.2[) and the Friedmann equation (|2.5|) take the form 0, 

d 2 d> . dd> ,. 



(2.10) 



h 2 (r) = \ 



(2.11) 



These coupled nonlinear differential equations completely define the time evolution of the inflaton field and the scale 
factor once the initial conditions are given at the initial time To. Namely, the initial conditions are fixed by giving 
two real numbers, the values of 4>{tq) and d(f>(To)/dr. 

It follows from eqs. (|2.11[) that 



d 2 a _ 1 
dr 1 ~ 3 



v{cj>) 



dr 



P+~ 3 P 



(2.12) 



When d 2 a/dr 2 > the expansion of the universe accelerates and it is then called inflationary. 
The derivative of the Hubble parameter is always negative: 

2 



dh 



1 

2 V dr 



(2.13) 



Therefore h(r) decreases monotonically with increasing r. Conversely, if we evolve the solution backwards in time 
from To, h(r) will generically increase without bounds. Namely, at some time r = r*, h(r) can exhibit a singularity 
where simultaneously a(r„) vanishes. 

In fact, the equations (I2.11j) admit the singular solution for r —> r*, 



(r) 



2 . T - T* 

3 log_ T- 



h(r) 



log a(r) 



1 



dr 3 (t — r* ) 

where & is an integration constant. The energy density e(r) and equation of state take the limiting form, 

1 P(r) 7 



P(r) 



3 (r - r*) 



-oo 



1 



(2.14) 



(2.15) 
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Namely, the limiting equation of state is p 
We have in this regime 



o(t) = C(t-t*)s ->0, (2.16) 



where C is some constant. That is, the geometry becomes singular for r — > t*. The behaviour near r* is non- 
inflationary, namely decelerated, since 

d 2 a r->r» 2 5 . . 

^2 = -g( T - r *) 3 ^-°°- ( 2 - 17 ) 

For r — > r», near the singularity, the potential u(^>) becomes negligible in eas. (|2.1ip . Therefore, eas. (|2.14p - (!2.17|) arc 
valid for all regular potentials v((p). 

The evolution starts thus by this decelerated fast-roll regime followed by an inflationary fast-roll regime and then 
by a slow- roll inflationary regime [2] . Recall that the slow-roll regime is an attractor t 4j , and therefore the inflaton 
always reaches a slow-roll inflationary regime for generic initial conditions. We display in fig. [T]the inflaton flow in 
phase space, namely d(f>/dr vs. <f> for different initial conditions. 

The number of efolds of slow-roll inflation N sr is determined by the time when the inflaton trajectory reaches the 
red quasi-horizontal line of slow-roll regime [see fig. [T]. We see that d<p/dr decreases steeply with <j>. This implies that 
N sr is mainly determined by the initial value of <fr with a mild (logarithmic) dependence on the initial value of difi/dr 

The inflaton flow described by eg. (12.141) results 

*< T > = V5 —i— < 2 - 18 > 

which well reproduce the almost vertical blue and green lines in fig. [TJ 

The inflationary regimes are characterized by the slow-roll parameters e v and r\ v [2] 

*-h{t) ■ 

The slow-roll behaviour is defined by the condition e v < 1/N. Typically, e v < 1/N during slow-roll. More generally 
accelerated expansion (inflation) happens for e v < 1 while we have decelerated expansion for e v > 1 as follows from 
eqs. (i2~TU|) - ([2~T^|) and pHty . 

The parameter r) v is also of the order 1/N during slow-roll and it is generically of order 1/N during fast-roll except 
when the potential v(<j)) vanishes. 

Eq. (|2.13p implies a monotonic decreasing of the expansion rate of the universe. There are four stages in the 
universe evolution described by eqs. (|2.11[) : 

• The non-inflationary fast-roll stage starting at the singularity r = r* and ending when d 2 a/dr 2 becomes positive 
[see eq. (|2~T2|) ]. 

• The inflationary fast- roll stage starts when d 2 a/dr 2 becomes positive and ends at r = T trans when e v becomes 
smaller than 1/N [see eq. ((2~T^) ]. 

• The inflationary slow-roll stage follows, and it continues as long as e v < 1/N and d 2 a/dr 2 > 0. It ends when 
d 2 a/dr 2 becomes negative at r = r en d- 

• A matter-dominated stage follows the inflationary era. 

The four stages described above correspond to the evolution for generic initial conditions or, equivalently, starting 
from the singular behaviour eqs. (|2.14j) . In addition, there exists a special (extreme) slow-roll solution starting at 
t = — oo where the fast-roll stages are absent. We derive this extreme slow-roll solution in sec. IIII A[ 

As shown in refs. 0, 0] the double well (broken symmetric) fourth order potential 
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FIG. 1: The complete inflaton flow in phase space, difi/dr vs. (f> for different initial conditions. We see that the inflaton always 
reaches a slow-roll regime for generic initial conditions represented by a red quasi-horizontal line. Hence, the slow-roll line is an 
attractor. Ultimately the inflaton reaches asymptotically the absolute minima dcfr/dr 
number of efolds of slow-roll inflation N ar increases for decreasing initial 
trajectories corresponding to N sr > 63 are colored in green. 



0, <j> = <p min = ^/8N/^ = 19-52 . . .. The 
> when dcfi/dr > initially. The (/> > 0, d(f>/dr > 



provides a very good fit for the CMB+LSS data, while at the same time being particularly simple, natural and stable 
in the Ginsburg-Landau sense. This is a new inflation model with the inflaton rolling from the vicinity of the local 
maxima of V((p) at <p — towards the absolute minimum <p = m/yX. 

The inflaton mass m and coupling A are naturally expressed in terms of the two relevant energy scales in this 
problem: the energy scale of inflation M and the Planck mass Mpi — 2.43534 10 18 GeV, 



M 2 
~Mp~i 



A = 



8iV 



M A 4 
Mpi) 



(2.21) 



Here N ~ 60 is the number of efolds since the cosmologically relevant modes exit the horizon till the end of inflation 
and y ~ 1 is the quartic coupling. 

The MCMC analysis of the CMB+LSS data combined with the theoretical input above yields the value y ~ 1.26 
for the coupling 0, 0] ■ V turns to be order one consistent with the Ginsburg-Landau formulation of the theory of 
inflation |2j . 

This model of new inflation yields as most probable values: n s ~ 0.964, r ~ 0.051 @, 0|- This value for r is within 
reach of forthcoming CMB observations. For y > 0.431946 . . . and in particular for the best fit value y ~ 1.26, the 
inflaton field exits the horizon in the negative concavity region V"(tp) < intrinsic to new inflation 0. We find for 
the best fit @,0], 



M = 0.543 x 10 16 GeV for the scale of inflation and 



1.21 x 10 13 GeV for the inflaton mass. 



(2.22) 
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We consider from now on the quartic broken symmetric potential eq. f|2 . 20[) which becomes using eq. (|2.8[) 



v(J>) = ~ 



— <A 4 H where a = - . 

4 V 4g y 87V 



(2.23) 



We have two arbitrary real coefficients characterizing the initial conditions. We can choose them as b and t* [see 
cq. (l2.f4j) ]. A total number of slow- roll inflation efolds N sr ~ 63 permits to explain the CMB quadrupole suppression 
[liSlal- Such requirement fixes the value of b for a given coupling y. 

We integrated numerically eqs. (|2.f f I) with eq. (|2.14p as initial conditions. We find that b = 4.745272 . . . f 0~ 5 yields 
63 efolds of inflation during the slow-roll era for y — 1 .26, the best fit to the CMB and LSS data. We find that b is a 
monotonically increasing function of the coupling y for fixed number of slow-roll efolds. At fixed coupling, b increases 
with the number of slow-roll efolds. 

We display in fig. [5] b as a function of y and the number of slow-roll inflation efolds N sr . 




500 



FIG. 2: Left panel: the coefficient b characterizing the initial conditions vs. the quartic coupling y for N sr = 63 efolds of 
slow- roll inflation. Right panel: b vs. N 3r for y — 1.26. The preferred values y — 1.26 and N ar are highlighted in both panels. 

For this value of y and 63 efolds of inflation during the slow-roll, fast-roll ends by r = r trans = 0.2487963 .... In 
figures[3j we depict log a(r), log/i(r), 0(t), log |0(r)|, log [AT e v (r)] and p(r)/p(r) vs. r till a short time after the 
end of inflation. We define the time T enc i when inflation ends by the condition a(r en d) — which gives (r en d — t*) = 
18.2547816.... 



Furthermore, we study in this paper the curvature and tensor fluctuations during the whole inflaton evolution in 
its three succesive regimes: non-inflationary fast-roll, inflationary fast-roll and inflationary slow-roll. 

The equation for the scalar curvature fluctuations take in conformed time n and dimensionless variables the form 



drf- 



S n (k;r)) = 



where dr\ = dr/a(r), 



1 d 2 



Wk(v) = - t~2 and z (^) 
z drj z 



i(rj) d<j) 



h(rj) dr 



(2.24) 



(2.25) 



10 




FIG. 3: Time evolution during during the three eras: non-inflationary fast-roll, inflationary fast-roll and slow-roll and beyond 
the end of inflation (MD era), loga(r), log/i(r), 0(r), log |0(t)|, log[iV £jj(t)] and p(t)/p(t) vs. r. a(r) grows monotonically 
reaching 63 efolds by the end of inflation. h(r) diverges for r — > t* = —0.8499574 . . . according to eo. (|2.14[ ) and decreases fast 
during fast-roll (r < Tirana = 0.2487963 . . .). Then, h(r) decreases slowly during slow-roll as discussed in sec. IIIIBI We depict 
h(r) for short times (0 < r — r» < 0.3) in fig. [5] 4>(t) diverges for r — > r» according to eq. (|2.14|) and decreases fast during 
fast-roll becoming very small during slow-roll. After the fast-roll stage where the infiaton field grows according to eq. (|2.14|) . 
4>(t) slowly rolls toward its absolute minimum at (j> en d = \/9>N/y — 19.52 . . .. log[iV tv(r)\ vs. r — t*. We have that e„(r„) = 3 
according to eqs. (|2.14[) and (|2.19|l . e„(r) decreases fast during fast-roll becoming of the order 1/N. We define the end of 
fast-roll (and beginning of slow-roll) by the condition N € v (t) = 1 which gives Ttrans — t* = 0.2487963 .... The equation of 
state p(t)/p(t) fastly decreases during fast-roll from the value p/p = +1 for r — > r* [see eq. H2.15|> ] passing through p/p = —1/3 
at the beginning of fast-roll inflation [see eo. (|2.12| )]. r = r s = r* + 0.0573, and reaching p/p = — 1 by the beginning of slow-roll. 
p/p vanishes again near the end of slow-roll inflation by r en d = t* + 18.698 .... 



In cosmic time r, cq. (|2.24l) takes the form 



d2 s d fc2 t r , ; 



S n (k;T)=0 



(2.26) 



where 



^(t) = 2 , : = ft (r) 
a 2 r 



2-7e v + 2e 2 v - y/%£ -g-tr -r] v (3- e v ) 



h 2 (r) 



h 2( \ r 7 , o 21 o ^ v '{<t>) „,,, 
h (r) [2 - 7 e„ + 2 e„J - 2 — T7 - v - u (0) , 



(2.27) 



dr h{r) 

and e„ and r\ v are given by eq. (|2.19p . 

We display Vr,(t) vs. t in fig. S]for the best fit value of the coupling y = 1.26 and 63 efolds of slow-roll inflation 
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The equation for the tensor fluctuations take in conformal time r\ and dimensionless variables the form Q 



S T (k;rj) + 



,2 a"{rj) 



a(rf) 



S T (k;r 1 )=0 



(2.28) 






-0.1 
-0.2 
-0.3 
-0.4 
-0.5 
-0.6 





/ / 




/ / 
/ / 


\ ( T - 

v\ 


1 / v rW ,'/ - 


\\ 


/ / 


\\ 


/ / 


\\ 


/ / 


\\ 


// 


\\ 


// 


\ \ 








\ \ 


-V 




/ 



16 18 



0.02 0.04 0.06 0.08 0.1 
X - X 




FIG. 4: The potential Vr.(t) felt by the fluctuations. Upper plot: Vjz(t) vs. (t — t*) in the stage where Vn(r) is repulsive 
{Vjz(t) > 0) which happens for (r — r*) > 0.114. Notice that Vizir) slowly decreases during the slow-roll stage as Vh(t) ~ 
2 h 2 {r) + 1 + 0(1/N) according to eq. (|2.27[ ) and fig. [3] Lower plots: Comparison of the exact (numerical) evolution and the 
analytic approximations eq. (|2.43[ ) during fast-roll and slow-roll. Left lower plot: (r — r») 2 Vtz{t) vs. r — r, in the stage where 
Vtz(t) is atractive (Vjz(t) < 0) from the exact (numerical^ calculation and from the analytic approximation eq. H2.43p . This 
happens for < (r — r») < 0.114. Notice that lim (r — t„) Vtz(t) — —1/9 according to eo. (|4.1[) . Lower right plot: Vtz(t) vs. 

T— >T* 

r — r* when Vtc(t) > from the exact (numerical) calculation and from the analytic approximation eg. (|2. 4311 . 



B. Inflaton and scale factor behaviour near the initial mathematical singularity 

In order to find the behaviour of </>(r) and a(r) near the initial singularity we write 

<^(r) = y|logl^4-0 1 (r) , h (r) = ^ \ r +h 1 (r). 



(2.29) 
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FIG. 5: Comparison of the exact (numerical) evolution (blue continuous line) and the analytic approximations (red dashed 
line) eq. (|2.40| l during fast- roll and slow-roll for log a(r), log h(r), 4>(r), logj0(r)|, e v (r) and p(r)/p(r) vs. r — t*. The exact 
In a(r) and \xih(r) are close to the approximation eq. (|2.40[ ). The scale factor is normalized to unit at r = 0, sixty efolds before 
the end of inflation. The exact (numerical) equation of state p(r)/p(r) is quite close to the analytic approximation eq. (|2.42[) 
both during fast-roll and slow-roll. The same happens for the exact (numerical) inflaton field 4>(t), 4>(t) and the analytic 
approximation eq. (l2.37|l . 



Inserting now eqs. (|2.29[) into eas. (|2.14p yields for 4>i( T ) an( i ^i( r ) the non-autonomous differential equations 

\ 3 



1 



3 h x 



V6 



hi 



2 T - T* 

3 



2 T T* 

3 







(2. 



3 (r — t*) 



6 



3 T — T* 



2 , r-r* 
log 

3 B 6 



J7_ 

12 



2 , T - T* 

3 



12 5 



where <fr stands for dip/dr. 

The asymptotic solution of eqs. (|2.30[) for t — > turns to have the dominant form 



h{r) T ^ T * (r -r,) 2 Pf (log- 



hi(r) 



(r - r*) Pi 1 lo. 



(2-31) 



where Pf(z) and P^(z) are fourth degree polynomials in their arguments. The polynomials turn to be of fourth 
degree because the inflaton potential is of fourth degree. Their explicit expressions follow after calculation 



V6 



( 10R' 4 

18 



1 . T — T* 7 



2, 3 t-t., 11, ,t-t, 49, r-r, 439 

- - log log 1 log 

3 8 6 3 6 b 9 B 6 54 

1 " 



13 



hi(r) 



log 



2 , r-r* 1 3 

- log - + — 

3 o 9 4 g 



log J 



log 2 



11 , t-t* 146 
— log 1 

3 8 b 9 



As a consequence, the scale factor near the singularity takes the form 



a(r) T = T * C (t - r*) 3 



1 + (r - t^) P 4 Q log 



(2.32) 



(2.33) 



where the coefficients of the fourth order polynomial P£ can be obtained from eqs. (|2.14p and ()2.32j) . 

C. Quantum loop effects and the validity of the classical infiaton picture 

When t — > quantum loop corrections are expected to become very large spoiling the classical description. More 
precisely, quantum loop corrections are of the order (H/Mpi) 2 Q. From eqs. (|2.7[) and (|2.14p the quantum loop 
corrections are of the order 



/ H \ 2 (t-t»)«1 

\Mp~l) 



3 (t-t.) Mpi 



1.66 10' 



1 / Tpi anc k 

9 \t-t, 



where we used m = 1.21 10 13 GeV 0. 

The characteristic time is here the Planck time 



Tpia 



m tpi 



M 



2.703 10~ 43 sec x m = 4.97 10" 6 . 



pi 



Namely, the quantum loop corrections are less than 1% for times 

. 10 

(r - r*J > — T Planck = 1.66 10 



(2.34) 



Therefore, for times (r — t*) > 10~ 5 the classical treatment of the infiaton and the space-time presented in sec. [TT] 
and III Bl can be trusted and we see that the classical description has a wide domain of validity. 

The use of a classical and homogeneous infiaton field is justified in the out of equilibrium field theory context as the 
quantum formation of a condensate during inflation. This condensate turns to obey the classical evolution equations 
of an homogeneous infiaton [ll[ . 

We see from eq. (|2.14p that the infiaton field becomes negative for t — > r*. But since a condensate field should be 
always positive, the classical and homogeneous infiaton picture requires 

T T* > b 

For the best fit coupling y = 1.26 and 63 efolds of inflation we have b — 4.745272 . . . 10~ 5 = 9.55 Tpi anc k which is 
consistent with eq. (|2.34[) . By comparing this value of b with eq. (|2.34[) we see that the quantum loop corrections are 
negligible in the stage where the condensate is already formed. 

We can obtain a lower bound on b since b increases with the number of inflation efolds N sr at fixed infiaton potential 
and since N sr cannot be smaller than the lower bound provided by flatness and entropy [2j. 

Although all inflationary solutions obtained evolving backwards in time from the slow-roll stage do reach a zero of 
the scale factor, such mathematical singularity is attained extrapolating the classical treatment where it is no more 
valid. In fact, one never reaches the singularity in the validity region of the classical treatment. In summary, the 
classical singularity at r = r* is not a real physical phenomenon here. 

The classical description with the homogeneous infiaton is very good for r— r* > 10 Tpi anc k well before the beginning 
of inflation. 
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D. The fast-roll regime: analytic approach 



As we see from fig. [3] the inflaton field 4>(t) is much smaller than dcfi/dr during fast-roll. We can therefore 
approximate the coupled inflaton evolution equation and Friedmann equation eqs. (l2.11[) as 



^ m *h d<t> n 
d^ +3h ^=° 



h 2 (r) = \ 



1 (d(j> 



2 V dr 



1 



(2.35) 



Or, in a compact form, 



which has the exact solution 



dr 



f0 
~d^ 



1 



3 dcj) 
2 dr" 



dr 



t — r* 



(r) 



Ti sinh 

where n turns out to be the characteristic time scale 

ri = 2 



log 



2ri 



tanh 



2ri 



(2.36) 



(2.37) 



6 TV 



We find for the best fit to CMB and LSS data, y = 1.26 and TV = 60, 

7i = 0.0592 = 11910 T P i anck , 



(2.38) 



(2.39) 



well after the Planck scale Tpi anck = 4.97 10~ 6 . 

The integration constant in eq. (|2.37j) matches with the small r — r* behaviour eq. (|2.14p . The Hubble parameter 
and the scale factor are here 



/i(r) = - — cothu , a(r) = C [n sinhw] 

3 Tl 



T* 



(2.40) 



where the integration constant was chosen to fulfil cq. (l2.16[) . The scale factor eq. (|2.40p interpolates between the 
non- inflationary power law behaviour eq. (|2.16p for t — t* — > and the eternal inflationary de Sitter behaviour for 
t — t„ 3> Ti. Since we have set v{<fi) equal to constant, slow-roll De Sitter inflation never stops in this approximation. 
Namely, neither matter-dominated nor radiation-dominated eras are reached in this approximation. 



We can eliminate the variable u between </> and d(j}/dr in eq. (|2.37[) with the result 



dr 



0(7 



(2.41) 



This equation generalizes eq. (|2.18|) which corresponds to the first term here and describes the behaviour for t — t* 
Notice that 



-co < 6(t) < 



2ri 



< — < +oo 

OT 



and that b/[2 Tl ] = 4.0105 10" 4 . 

The evolution described by eqs. (|2.37|l - (|2.40p starts from the mathematical singularity at t = r» with monotonically 
decreasing dcfy/dr and h(r) and a monotonically increasing 4>(t) from its initial value 4>(t*) — — oo. 

Slow-roll is reached asymptotically for large t since d(f>/dr vanishes for t — t„ — > oo. 
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We find for the parameter e v [eq. (12.19[) ] and for the equation of state, 

3 p(r) 2 



e„(r) 



1 + sinh u 



P(r) 



cosh 2 ■ 



(2.42) 



We see that e v (r) monotonically decreases with r and vanishes for r — T. t — > oo. The equation of state p/p smoothly 
interpolates between +1 at t = t, (extreme non-inflationary fast-roll) and —1 (slow-roll inflation) for r — t* — > oo, 
passing by p/p = —1/3 (the beginning of fast-roll inflation) at r — r* = 0.0573. 

The potential Vr(t) eq. (|2.24[) felt by the fluctuations takes here the form 



Vn{r) = ±- 
65 



1 



1 



2 sinh u cosh u 



Tl 



(2.43) 



The limiting values of h{r), 4>(t) and Vh{t) for r — > 00 give a reasonable approximation to the numerical results. 
We have 



h(oo) 




(j){oo) 



log 



2ri 



(00) = , V w (oo) 



4 AT 

^7 



(2.44) 



The characteristic time scale T\ is generically a small number since according to eq. (|2.38p T\ ~ 1/yN. The value of 
ti for the best fit value for y is given in eq. (|2.39l) . 

The end of fast-roll Ttrans can be estimated in this approximation by using eq. (|2.42|) for c v (t) setting e v (T tr ans) = 
l/N. This gives, 



e v (t) ~ 12 e" 



1 



1 



ti ln(12 AT) = 0.195 . 



0.2487963 . 



This approximated value for T tra ns should be compared with the exact numerical result T trani 
h(jt ra ns) and Vn(Tt ra ns) differ in less than 1% from their values at r = 00 given by eq. flZBJ) . 

In figs. [5]we plot lna(T), ln/i(r), </>(t), In |<^(t)|, c v (t) and p{t)/p{t) computed numerically and computed 
using the analytic expressions eqs. (|2.37|) - ()2.42j) . We compare in figs. [4] the exact potential Vr(t) with the analytic 
approximation eq. (|2.43p . 

We see that the simple analytic formulas eqs. (l2.37l) - (|2.43l) provide a very good approximation during the fast-roll 
regime r < ttrans = 0.2487963. . .. In particular, eq. (|2.37|) provides an excellent approximation to </>(t) as shown in 
fig. [5] In particular, the analytic formulas eos. (|2.37[) - (|2.43p become exact near the singularity at r = t*. 



E. The fast-roll regime: numerical solution 

To construct a singular solution we can integrate eqs. (|2.11j) backwards in time starting from initial conditions of 
strong non-inflationary fast-roll type, namely 

producing a given total number N sr of slow-roll inflationary cfolds. For instance, we start from some <f> and <j> such 
that K = 10 4 . The time extent backwards from this moment has to be limited so that, integrating back and forth, 
the required relative accuracy of 10~ 12 is preserved. We furthermore impose that A^ sr = 63. 

We adopt the convention that conformal time 77 vanishes from below when inflation ends and that <z(t = 0) = 1 
when there are still A^ = 60 efolds till the end of inflation. This choice of the scale factor normalization seems the most 
natural. Then, r\ has a finite non-zero limit 77* as r approaches the time r* of the singularity, since air) ~ C (r — t*) 1 / 3 
as t — » t* according to eq. (|2.16[) . That is, 

dr' r T dr 1 



1= ^ =1] * 



l Tend a(r>) " J Tt a(r') " 
The numerics of a fast-roll solution of this type are in Table I where a relative accuracy of 10 -12 is preserved. 
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Using the asymptotic behaviour eq. (12.14[) as r — !• r+ we obtain from Table I: 

r* = -0.8499574 .. . , b = 4.745272 ... 1(T 5 and 77, = -15.605614 ... . 

Slow-roll begins at T trans = r* + 0.2487963 . . . = -0.6011611 . . .. 
The initial value of the ratio 

dip/dt <f> 



has the dimension of mass. The natural mass scale in the problem is here the energy scale of inflation M. Therefore, 
assuming this ratio of the order M yields 

Z < 10 3 . 

<p m 

Hence, it s natural to start the fast-roll evolution with <j)/(f> < 10 3 . 





K = 5.3458 ... 10 7 


K = 10 4 


inflation start: a = 


fast-roll — > slow-roll 


a = 1 


inflation end: a — + 


T 


-0.8499493 . . . 


-0.8493593 . . . 


-0.7746494 . . . 


-0.6011611... 





17.4048242 . . . 


<t> 


-1.4401237. . . 


2.0690604 . . . 


5.9342489 . . . 


6.4783577 . . . 


6.7484076 . . . 


18.5586530 . . . 


<j> 


100391.035 . . . 


1365.05241 . . . 


8.8601670 . . . 


0.9182661 . . . 


0.3974015 . . . 


0.94150557... 


log a 


-7.0325621 . . . 


-5.5999353 . . . 


-3.9142151 . . . 


-2.9999999 . . . 





60 


h 


40984.4689 . . . 


557.30817... 


6.2650841 . . . 


5.0295509 . . . 


4.9653990 . . . 


0.6657449 . . . 


V 


-15.6050091... 


-15.376218... 


-15.3549996... 


-4.0169827 . . . 


-0.2020609 . . . 






TABLE I: Fast-roll solution with jV sr = 63 efolds of slow-roll inflation. Recall that r = 4.97 10" 6 (t/t P i anck ). 



III. THE SLOW-ROLL INFLATIONARY ERA 



The extreme slow-roll solution 



There always exist a special solution of eas. (|2.1ip that starts at r = —00 with vanishing inflaton, vanishing scale 
factor but nonzero Hubble parameter. More precisely, eqs. (l2.1ip can be approximated for small <j> and <p as 



dr 2 dr 
h 2 (r) - I v(0) 



0-0 , 



(3.1) 



where we used eqs. (|2.9l) and (12. lip . 

Eqs. p.ip admit the asymptotic solution for r — > — 00 



, e 



, %) = 



v(0) 



o(r) 



-00 

: e V 3 



where Cq is an integration constant, v(0) — 2N/y for the double-well potential eq. (|2.23p and 

1 



y/3 v(0) + 4 - y/3 V (0) 



> 



(3.2) 



Notice that a can be expressed in terms of the fast-roll characteristic time-scale t\ [eq. (|2.38l) ]. 
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since n ~ 0.0592 < 1 [see ea. (j239)) ]. 

It must be noticed that the characteristic time scale of the inflaton evolution in the extreme slow-roll solution for 
early times [see eq. (|3.2[) ] 

a ri 

turns to be the inverse of the characteristic time scale t\ of the fast-roll solution and to be very large. 
On the contrary, the characteristic time scale of the scale factor evolution in the same regime is very short 



«(0) 



= 3 n -c l 



The fast-roll stages both non-inflationary and inflationary are absent in this solution. The extreme slow-roll solution 
only possesses the slow-roll inflationary stage followed by the matter dominated era. 





inflation start 


a = 1 


inflation end: a = + 


T 


-344.9514017 . . . 





17.40482446 . . . 




1(T 8 


6.7484118 . . . 


18.5586530 . . . 


4> 


alO" 8 = 5. 8937108453. ..10" 10 


0.3973384 . . . 


0.94150557... 


log a 


-1938.4867948... 





60 


h 


y/2N/(3y) = 5.6361006... 


4.9653973 . . . 


0.6657449 . . . 


V 


— oo (f.a.p.p) 


-0.2020610 . . . 






TABLE II: Relevant quantities of the extreme slow-roll inflaton solution for the coupling y = 1.2592226.... We adopt the 
convention that a(r = 0) = 1 when there are still N = 60 efolds till the end of inflation. Recall that r = 4.97 10 -6 (t/tpianck)- 



For the value of the coupling y = 1.2592226 . . ., we get for the extreme slow-roll solution 
a = 0.058937108 . . . , </> ond = 18.5586530 . . . , <p cnd = 0.9415055 . . . 
In table II we display the values of the relevant magnitudes for this extreme slow-roll solution. 



(3.3) 



Except for the extreme slow-roll solution, all solutions are of fast-roll type and come from singular values of <fi and 
h according to eq. (|2.14p as r —> r+ for some finite t* characteristic of each particular solution. The slow-roll stage 
(which starts when e v = 1/N from above, and ends when again e v = 1/N from below) of all distinct solutions turns 
to be almost identical to that of the extreme slow-roll case as one could expect for an attractor. 



B. The inflaton during slow-roll inflation: analytical solution 

In the slow- roll regime higher time derivatives can be neglected in the evolution eqs. (|2.1ip with the result 

3*(r)0 + t/(0) = O , h 2 (T) = ^±. 
These first order equations can be solved in closed form as 



N[<f>] = - I v(<f>') ^- 



(3.4) 



(3.5) 



where iV[</>] is the number of e- folds since the field <fi exits the horizon till the end of inflation (where it takes the value 

(t>end)- 

Eq. (|3.5j) indicates that N[(j)] scales as (ft 2 and hence the field 4> is of the order y/N ~ y/60- Therefore, we proposed 
as universal form for the inflaton potential [H, [l(| 



v{<p) = N M 4 w( X ) , 



(3.6) 
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where x is the dimensionless, slowly varying field 

X = 

The equations of motion (|2.11[) in the field x become 



N M 



pi 



N 



(3.7) 



UHt) 



1 d 2 x 



1_ (dx 

2 N \ df 



^? hz + w (x) 



with % 



iV 



(3.8) 



and f stands for the rescalcd dimensionless time 



m t 



N VN 

To leading order in the slow-roll approximation (neglecting 1/7V corrections), eqs. p.8[) are solvable in terms of quadra- 
tures 



T - T tr 



where Ttrans stands for the beginning of slow-roll inflation and we used that 



HI 



w(x) 



For the broken symmetric potential eq. (|2.20p . from eqs. (|2.10[) , (|3.9[) and (|3.10p . we find 

X(f ) = Xinrans) + O (I) = e"V¥ + O (i 



1 - e ~V^? (fc " d ~ 

1/ 1 



W(f) = 

-n = -1 + 

p lTJ 6 N sinh 2 [7f (f end - f) 



O 



1 

JV 



1 



for f irona < f < f end = \ — In 

2 y 



8 



X 2 (ftrans) J/ 



(3.9) 

(3.10) 

(3.11) 

(3.12) 
(3.13) 



Inflation ends when the equation of state becomes p/p = —1/3 [see eq. (|2.12j) ]. According to eq. (|3.12j) . this happens 
when T en d — f ~ O (\/^JTt*j . Therefore, expressions eqs. (|3.11|) - (l3.12|) are valid as long as 



Ttrans ^ T ^ Tend & 



where O 



> . 



That is, eas. p.lll) hold while the inflaton is not very near the minimum of the potential Xend — \/8/y. 
By integrating the Hubble parameter H{t) we obtain for the scale factor a(f ) 



off) [~T Ar ., „ . TV 2/A . 

log —~ r = \ — N (T - Tfrans) ~ — X (Ttrans) 

a(Ttrans) \ oy 8 



(3.14) 



\2N 1 

— Tfl \ t ttrans ) 77 

3 V 8 



^(^trans ) 


2 


Mpi 





e V?l- m (*-*«ron.) _ ^ 



where we used eqs. (|2.7l ) and (|3.1ip . It must be noticed that a(f) is not exactly a de Sitter scale factor, even in the 
large N limit at fixed f. 



19 



At the end of inflation the number of efolds is In a ~ 64, the infLaton is near its minimum 

■ 2.52 , 

X starts to oscillate around zero and H(f) begins a rapid decrease (see figs. [3j. At this time the inflaton field is no 
longer slowly coasting in the w"(x) < region but rapidly approaching its equilibrium minimum. When inflation 
ends, the inflaton is at its minimum value up to corrections of order 1/%/iV. Therefore, we see from the Friedmann 
eq. iprgjt and eqs. (f3~TTj) that 

Jf (^) ^end) = O (J^j , w( X (f end )) = O (Jj) and therefore, H(f end ) = O (^=j , (3.15) 



while H(f trans) = 0{\). Namely, the Hubble parameter decreases by a factor of the order \JN ~ 8 during slow-roll 
inflation. We see in fig. [3] that the exact H(t) decreases by a factor six during slow- roll inflation, confirming the 
slow-roll analytic estimate. 

We can compute the total number of inflation efolds N to t to leading order in slow-roll inserting the analytic formula 
for f end eq. (|3.13[) in eq. (|3.14p with the result, 



Ntot = — < In 

y 



X 2 (f trans) V 



l + \yX 2 (nrans)\+0(-^= ) . Ci.K,) 



We have verified the slow-roll analytical results eqs. (|3.1ip - (|3.16p comparing them with the numerical solution of 
eas. (|2.11j) . Both results are concordant up to the error estimation in each case: O (1/iV) or O ^1/Va"^ 

The field <j> as a function of the dimensionless time r eq. (|3.1ip takes the form 

4>{T) = 4>{Ttrans) e V 7 ^ ) 

and then 



For y ~ 1.26 and N = 60 we get \/y/[6 N] = 0.0577 in agreement with the slope of the red quasi-horizontal slow-roll 
line in the phase space flow fig. [TJ 

IV. COMPLETE FLUCTUATIONS EVOLUTION AND FAST-ROLL EFFECTS ON THE POWER 

SPECTRUM. 

A. Scalar and tensor fluctuations near the initial singularity. 

In order to study the curvature and tensor fluctuations in this regime, it is important to evaluate the parameter e v 
and the potential felt by the fluctuations Vr. 

Inserting eqs. (|2.29p and (|2.3ip into eas. (|2.19p and (I2.27P yields near the initial singularity 



V n (r) T ^= ' 



9(r-r,) 2 



1 + (r - t*Y Pi log 



e v = 6 



1 + {T-T*y PI lO; 



(4.1) 



Wniv) V =° ~ [1 + ^ 3 Pf (log V)} ■ (4-2) 



where 



T^T„ 3 , .2 

is the conformal time for r — > r* and P^ (x), P£(x) and P^ (x) are polynomials of degree four in x 
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We see that the fluctuations feel a singular attractive potential near the r\ = singularity. Actually, the behaviour 
of Wn(r)) for 77 — ^ is exactly the critical strength (—1/4) for which the fall to the centre becomes possible in a 
central and attractive singular potential 

We find from eas. (|2.26[) and (|4.1j) for the fluctuations near the singularity 



An{k)+B n {k) log J- 
Vo. 



(4.3) 



where rjo is the time when the initial conditions will be imposed, A-jz and B-n are complex constants constrained by 
the Wronskian condition (that ensures the canonical commutation relations) @ 



W[S n ,S* n ] = S n ^-^S 



drj dr\ 



■r 



Namely, 



2 Im[A n B* n ] 



(4.4) 



(4.5) 



Precisely, the logarithmic behaviour for 77 — > of the wave function eq. (|4.3j) describes the fall to r — r* =0 for 
the critical strength of the potential W-R.{rj). For larger attractive strengths the wave function ea. (|4.3l) shows up an 
oscillatory behaviour Q. Notice, however the physical nature of the process: here we have a time evolution near a 
classical singularity at a given time while in the potential case one has particles falling (or emerging) from a point in 
space where the potential is singular. 

In general, the mode functions for large k must behave as free modes (plane waves) since the potential Wji{rj) 
in eq. (|2.24[) becomes negligible in this limit except at the singularity r = r*. One can then impose Bunch-Davies 
conditions for large k which corresponds to assume an initial quantum vacuum Fock state, empty of curvature 
excitations [2| 



Sn(k;r) 



-i k 77 



<2 k 



(4.6) 



and therefore 



^j^(fc;r?o) fe =°° -i k S n (k;r] ) . 
drj 

Eq. (l4.6[) fulfils the Wronskian normalization eg. (14.41) . 

In asymptotically flat (or conformally flat) regions of the space-time the potential felt by the fluctuations vanish 
and the fluctuations exhibit a plane wave behaviour for all k (not necesarily large). This is not the case near strong 
gravity fields or curvature singularities as in the present cosmological space-time where W-r (77) can never be neglected 
at fixed k. However, we can choose Bunch-Davies initial conditions (BDic) at 77 = 770 by imposing 



drj 



(k; rjo) = —i k Sn(k; rjo) for all k . 



(4.7) 



That is, we consider the initial value problem for the mode functions giving the values of S-R,(k;r)) and dSu/dr] at 
V = Vo- 

Notice that eq. (|4.7j) combined with the Wronskian condition eq. (|4.4[) implies that 



\Sn(k;r) )\ 



1 



f2k 



which is equivalent to eq. 



for large k. 



dS n 



drj 



(k;m] 



Since the mode functions Sn{k\rj) are defined up to an arbitrary constant phase we can write eq. (|4.3l) valid near 
the metric singularity as 



?7,77o->-0 



1) 



2 k i] Q 



I ( l+ikr)o) log^- 
2 J Vo 



(4.8) 
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Inflation (fast-roll) starts 


a(r a ) = 


t s = r, + 0.0753090 


Vtz(t) becomes positive 


Vk(t+) = 


r+ = t» +0.114 


End of fast-roll 


N e v (Ttrana) = 1 


n„„ s = r, + 0.2487963 . . . 


Maximum of Vjz(t) 


V£(r M ) = 


tm = t» + 0.3503 , V^(tm) = 51.196 


End of Inflation 


p(r e nd) = 


r e „ d = t» + 18.2547816 



TABLE III: Selected time values for y — 1.26 and AT sr = 63 efolds of slow-roll inflation. Notice that slow-roll starts exactly 
when fast-roll ends. Recall that r = 4.97 10~ 6 (t/tpi anc k). 



In eq. (|4.3j) this corresponds to the coefficients, 

1 



An = 



/ 2k 



B 



- + i k r) 



Tk \2 



We have in cosmic time, 



'2 k 



to 



1 - - + i k T 3 log 



TQ 



(4.9) 



where t — r* = (2 rjo /3) 5 for r — > r* . 

Namely, imposing the BD initial condition (BDic) eq. (|4.7p at small 770 where the small r\ behavior eq. (|4.3[) applies, 
yields specific values for the coefficients of the linearly independent solutions y/rj and ^Jfj log 77 that we can read from 
eas.(01l-(01l. 



For general tq (i. e., tq not near r*), the mode functions for t — > r* take the form 



Sn(k]T) 



'2 k 



TO 



X(k,T )- Y(k,T ) + 



i k 



X(k,T ) 



log- 



(4-10) 



where we imposed eq. (l4.5[) and we have from eq 



X(k,T*) = l and Y(k,T*) = ~ . 

Notice that X(k, to) > for tq — > r* as we see from eq. (|4.9|) . Our numerical calculations show that X(k, to) > 
for all t and k. 



B. The primordial power spectrum, Scalar curvature fluctuations and the CMB+LSS data. 



The power spectrum of curvature perturbations 1Z is given by the expectation value < 1Z 2 > in the state with 
general initial conditions [2] 



< K 2 (x,r)) >-- 



( m 



PL 



\S n (k;v)\ 2 k 2 dk 



z 2 (r)) 2tt 2 



(4.11) 



where z(t]) is given by eq. ()2.25|) . Notice in eq. (|4.1ip the factor (m/Mpx) 2 in the physical power spectrum expressed 
in terms of the dimensionless quantities used here. 

The power spectrum at time 77 is customary defined as the power per unit logarithmic interval in k 



f°° dk 

<n 2 (x, V ) >= — P n {k,v) 



Therefore, the scalar power for general initial conditions is given by the fluctuations behavior by the end of inflation 



1, 



Pn(k) = 



( m 



2 fc 3 



i 1 , li m 

V MpL, J 2 TT 2 r^Q- 



Sn{k;v) 



z{rj) 



(4.12) 
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The mode functions S-ji(k;rj) obey the fluctuations equation (|2.24[) where the potential W-n{rf) [eg. (12.271) ] during 
slow-roll and to leading order in 1/N takes the simple form [2j , 



WnW = - 

if 



1 + 2 ( 3 e « - ^) 



TC 4 



• = - + '■-'>< - + C° ( 



(4.13) 



In the slow-roll regime we can consider e v and r\ v [see ea. (|2.19|i ] constants in time in eq. (|4.13p . During slow- roll, the 
general solution of eq. (|2.24p is then given by 



with 



Sn(k; 1]) = A n (k) g u ^ (fc; rf) + B n (k) g* (k; ry) , 



(4.14) 



(4.15) 



Afi(k), B-ji(k) are constants determined by the initial conditions and Hv (z) is a Hankel function. 
The Wronskian of the solutions S-r, is given by eq. (|4.4l) and 

W[g„,g*] = i 

This generically determines that 

\A n (k)\ 2 - \B n (k)\ 2 = I . (4.16) 
For wavevectors deep inside the Hubble radius \kr}\ 1 the mode functions g l/ (k;r]) have the asymptotic behavior 

1 _„* i. ^ ^ /, \ n— >— oo 1 



while for r\ — > 0~, they behave as: 



77— > — oo 



-i k r) 



'2 k 



t)->o- r(i/) 



V2vr k \ikf] 



(4.17) 



(4.18) 



In particular, in the scale invariant case v = | which is the leading order in the slow-roll expansion, the mode functions 
eqs. (|4.15[) simplify to 



-i k rf 



2k 



1 



k rj 



As we see from eq. (|2.25|) . z[if) obeys eq. (l2.24[) for k — and therefore z(rj) in the slow-roll regime behaves as 

zo 



z{rf) = 



{-ho ri) u -*—i 



(4.19) 



(4.20) 



where zq is the value of z(r]) when the pivot scale fco exits the horizon, that is at i] = — 1/fco- Combining this result 
with the small r\ limit eq. (14.18|) we find from eas. (|4.12|) and (|4.20j) . 



Pn(k)=Pi D (k)[l + D(k)} , 

where we introduced the transfer function for the initial conditions of curvature perturbations: 

D{k) = 2 \B n {k)\ 2 - 2 Re [A K (k) B* n {k) i 2 ^} . 

D(k) is obtained imposing BDic at r = tq according to eq. (|4.7p . 

Notice as shown in sec. IV Al that the transfer function D(k) enjoys the properties 



(4.21) 
(4.22) 



1 + D(k) k =° 0{k n ° +1 ) , D{k) k ^° 



° [ k 2 



(4.23) 
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D(k) accounts for the effect in the power spectrum both of the initial conditions and of the fluctuations evolution 
during fast-roll (before slow-roll). D(k) depends on the time r at which BDic are imposed. 

If one chooses the extreme slow-roll solution presented in sec. IIII Al and imposes BDic at tq = — oo (that is, 
r/o = — oo) then D(k) = and the fluctuation power spectrum at the end of inflation is the usual power spectrum 

P n (k)=P£ D (k). 

P-jz D (k) is given by its customary slow-roll expression, 



log P£ D (k) = \ogA s (k ) + {n s 



^ , k 
1 log - 

fc 



fc 



(4.24) 



We solved numerically the fluctuations equation (|2.26j) in cosmic time with the BDic eq. (|4.7|) covering both the fast- 
roll and slow-roll regimes. We started at initial times tq ranging from the vicinity of t = till the transition time 
Ttrans = 0.2487963 . . . from fast-roll to slow-roll. We computed the transfer function D{k) from the mode functions 
behaviour deep during slow-roll inflation from eqs. (|4.12|) and (|4.21j) 2]. In figs. |6]we depict 1 + D(k) vs. k for twelve 
values of the time r where BDic are imposed. 
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0.004 
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0.007 
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Ax = 0.02 
Ax = 0.04 
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At = 0.08 



At = 


0.1 


At = 


0.15 


At = 


0.2 


At = 


0.2487 




1.5 1.5 

log(k/m) 



-1.5 1.5 

log(k/m) 



FIG. 6: Numerical transfer function 1 + D(k). Lower left panel: Numerical transfer function 1 + D(k) for BDic at r = to = 
t» + At, for different At values as given in the picture. We see here that the peak of 1 + D(k) grows and moves for larger k 
as to increases. Here N sr — 63. Lower right panel: The transfer function 1 + D(k) when the BDic eq. (|4.7[) are imposed during 
slow-roll at finite times to and N ar efolds of slow-roll have still to occur. Upper panels: Numerical transfer function 1 + D(k) 
for BDic at r = to = t„ + At, for different values of At as given in the picture. We get stronger oscillations in 1 + D(k) for 
decreasing tq in the range At < 0.04. Here N sr = 63. 



Notice that when BDic are imposed at finite times r , the spectrum is not the usual P-£ D {k) but it gets modified 
by a non-zero transfer function D(k) eq. (|4.2I p . The power spectrum Pji(k) vanishes at k — and exhibits oscillations 
which vanish at large k [see figs. [5] and [7]. 



24 



140 




BDatN = 


61 


sr 




BDatN = 


63 


sr 




BDatN = 
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log(k/m) 



FIG. 7: Power spectrum with BDic eg. (|4.7fl imposed during slow-roll when N ar efolds of slow-roll inflation have still to occur. 
We see here the decrease of the power spectrum Ptz as k" s ~ multiplied by the oscillations of 1 + D(k). See eas. (|4.21[l and 
(|4.24[) and figs. [6] The non-oscillatory black curve corresponds to the usual power with BDic at 770 = -co eq. (|4.24l ) decreasing 
as fc™ a_1 . The later are imposed the BDic, the smaller is the number of slow-roll efolds N sr and the whole fe-spectrum shifts 
to larger k. 



During slow-roll different initial times To lead essentially to a rescaling of k in D(k) by a factor 770 since the conformal 
time r\ is almost proportional to 1/0(77) during slow-roll [see figs. [TJ-H] and below eq. (|5.7l) ]. By virtue of the dynamical 
attractor character of slow-roll, the power spectrum when the BDic are imposed at a finite time To cannot really 
distinguish between the extreme slow-roll solution (for which slow-roll starts from the very beginning 770 = —00) or 
any other solution which is attracted to slow-roll well before the time To- 

C. Accurate numerical computation of the power spectrum and the transfer function D(k) of initial 

conditions. 



In order to accurately calculate n s we proceed as follow. We match the solution Sn(k; 77) with the slow-roll solution 
g VTi {k]rf) eq. (|4.15p at the time To when N sr efolds of slow-roll have still to occur. 77 and v-r. are computed at this 
time To. In practice, this corresponds to setting An(k) — 1, B-ji(k) — (and therefore D-ji(k) = 0) in the Bogoliubov 
transformation eq. (|4.14p . 

Then, we integrate numerically the fluctuations equations eq. (|2.26[) . By construction, this produces the standard 
spectra P^ D (k) eq. (|4.24[) that quickly stabilize as N sr is increased a few efolds above N = 60. 



It is convenient to introduce the quantity 



L„ = loe 



M PL 



A s (ko = 777) 



(4.25) 



with fco 
of N, r . 



m when 0(77) = 1, that is N — 60 efolds before inflation ends. In table IV we provide L s for several values 
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Nsr 


L 3 


n a 


Tlrun 


61 


4.6585381 . . . 


0.9637013 . . . 


-0.0000701 . . . 


63 


4.6583004 . . . 


0.9641135... 


-0.0001639 . . . 


65 


4.6584371 . . . 


0.9642483 . . . 


-0.0002165 . . . 


67 


4.6584463 . . . 


0.9642444 . . . 


-0.0002165 . . . 


69 


4.6584469 . . . 


0.9642448 . . . 


-0.0002167... 



TABLE IV: Exact values of L 3 = \og[{MpL /m) 2 A s (ko = m)] for several values of N ar from the numerical calculation. The 
exact values of n s vary little with N sr and are close to the slow-roll approximation value. Also n run is close to the value in the 
slow-roll approximation. 



To transform this fco m a wavenumber today we need: 

• the total redshift from 60 efolds before inflation ends till today [since we choose a(r = 0) = 1 when there are 
still N = 60 efolds till the end of inflation] . 

• the value of to as determined by the observed value of the amplitude A s (fco). 

Let kQ MC be the value of the pivot scale of CosmoMC [that is 50 (Gpc) -1 today] 60 efolds before the end of 
inflation. Then, we have from eqs. (|4.24p and (|4.25|) . 



log A s (k Q = to) = L s + 2 log -f- = L™ c + (n? MC - 1) log 



M PL ~ s ^ ^ fcCMC 



log 



{.CMC 



' 1 



where L^ MC = log A^ MC (kQ MC ) and n™ c are best fit values in a given CosmoMC run. Since the running index 
n run is O (l/N 2 ), we get for to, 



CMC_ 

1 



1 + 0, A- 



(4.27) 



The wavevectors at a = 1 (60 efolds before inflation ends) and today are related by [2| 

60 

k a=1 = — k today , (4.28) 
a r 

where a r is the scale factor by the end of inflation 



a r = 2.5 10- 29 , (4.29) 

V -"60 

and Hqq is the Hubble parameter 60 efolds before inflation ends. We thus have for the pivot wavenumber at a = 1 



and 



fc CMC ~ 1.46 . . . x 10 15 GeV (4.30) 



16.67... \ 8 , r CMC t \ v. i~ ^60 

exp [Lg — L s j , where h eo 



Notice the small 1/N correction (n^ MC — l)/2 in the exponent of m/Mp^. Eq. (|4.26l) yields for the best fit CosmoMC 
run i^MC = _ i9.9 8 08 . . . and n s = 0.9635 . . . @: 

to ~ 4.8114 . . . 10~ 6 M PL = 1.1717 . . . 10 13 GeV 

The exact values given above in Table IV 



-U = ( ™ \ exp(L s ) , n s and n r 
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are obtained taking into account the fast-roll and slow-roll stages in the numerical calculation. We can compare them 
to their slow-roll (leading 1/N) analytic counterparts for the double- well quadratic plus quartic potential, Q 

N* f m \ 2 (!-*)« y 3z+l y 2 z 

As -\2^\M PL ) y 2 z ' Us ~ 1 AT (1 — z) 2 ' Urun ~ N 2 (1 - z)^ 24 z 35 Z + 3 > 

where N — 60, z = 0.117446 and y = z - 1 - log z = 1.2592226 . . ., that is 
N 2 ( m \ 2 

A s = ——r — — exp(4.59536898 . . .) , n s = 0.9635620 . . . , n run = -0.0000664 . . . 
\2-k z \MpLj 

The figure in the exponent is to be compared with the L s values in Table IV. The agreement with Table IV is quite 
good, especially for n s . 




log(k/m) log(k/m) 

FIG. 8: Upper Left panel: Difference between the (approximate) transfer function D(k r/o) eqs. ()5.7[) - (|5.9l) for v-n — 1.5182189 . . . 
and the numerical (exact at least to a 10 -7 relative error) transfer function D(k), when N sr — 63. Upper Right panel: Difference 
between the (approximate) transfer function D(k yo) eqs. H5.7p - (|5.9p for v-r, = 3/2 (the scale-invariant value) and the numerical 
(exact) transfer function. We see that the difference in the right panel [eq. (|5.9|) ] is < 0.014 while in the left panel the difference 
of the analytic formula eq. (|5.7p is much smaller, < 0.0005. Lower Left panel: difference between the exact (numerical) D(k) 
computed for the fast-roll inflaton solution of table II and for the extreme slow-roll inflaton solution of table I when BDic are 
imposed 63 efolds before the end of inflation. Lower Right panel: difference between the numerical (exact) fast-roll D(k) and 
the approximate D(k yo) calculated with v-r = 3/2 and 770 = —4.0169827 . . .. We see that the differences are small in both 
cases. 
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We now find the exact (numerical) transfer function D(k) for the initial conditions, by simply taking in ea. (|4.21|) 
the ratio of the two power spectra: P-ji(k) with BDic at time To and P^ D (k). In the case of BDic at finite times the 
result is given in fig [6] At the largest value k/m = 100 of the wavenumber interval considered, we have 

1 + I>(100 to) = 0.9996994 . . . , 1.0000061 . . . , 1.0000001 . . . 

for 

N sr = 61, 63, 65 and 67, respectively. 

This provides a good check of the accuracy of the calculation. 

In figs. [S]we compare the numerically computed D{k) against D(k rjo) analytically computed for BDic imposed at 
time rjo during slow-roll in eq. (|5.7l) . sec. IV Bl The comparison is performed for BDic imposed when N sr = 63 on the 
extreme slow roll solution, which corresponds to i]o = —4.0202308 . . .. We consider two values of v-n : v-ji = 2 — n s /2 = 
1.5182189 . . ., n s — 0.9635620 . . . corresponding to slow-roll at leading 1/N order, and the exactly scale-invariant case 
v n = 3/2. Notice that in the latter case D(k 770) has the explicit simple analytic form eq. (|5.9l) . 

The maximum of the numerical transfer function 1 + D{k) is located at k/m = 0.68755... and has the value 
1.13218 . . . The maximum of 1 + D(k 770), when v-n — 3/2 is in k/m — 0.68755 . . . and has the value 1.13009 . . .. Recall 
that these values of k/m have the scale fixed by the choice a — 1 when N — 60 efolds lack before inflation ends. 

Let us now consider the fluctuations on the fast-roll solution of Table II. Since 77 has a finite lower limit, the choice 
Aiz(k) = 1, Bfi(k) — has little meaning and BDic can be imposed only at a finite time To later than the singularity 
time r*. If t is exactly the transition time T trans when e v = 1/N, fast-roll ends and slow-roll begins, (to proceed 
for N sr = 63 efolds), then D(k) does not differ too much from that computed with the extreme slow roll solution. 
This comparison is performed in the lower left panel of fig. [8] In the right panel D(k) is compared to the D(k rjo) for 
v n = 3/2 and 770 = —4.0169827 . . ., which is the value of the conformal time at the onset of slow-roll (see Table II). 

When the BDic are imposed during the fast-roll stage well before it ends, D(k) changes much more significantly 
than along the extreme slow roll solution. This is due to two main effects: the potential felt by the fluctuations is 
attractive during fast-roll and 770, far from being almost proportional to 1/0(77), tend to the constant value 77* as 
t — > r + and 12(77) —> 0. The numerical transfer functions 1 + D(k) obtained from eqs. (|4.12[) and (14.211) are plotted in 
figs. E] 

The fact that choosing BDic leads to a primordial power and its respective CMB multipoles which correctly 
reproduce the observed spectrum justifies the use of BDic for the scalar curvature fluctuations. 



D. The effect of the fast-roll stage on the low multipoles of the CMB 

In the region of the Sachs- Wolfe plateau for I < 30, the matter-radiation transfer function can be set equal to unity 
and the CMB multipole coefficients C[s are given by @ 



A f OO 11 

c i = YJ T Px(k) {mm ~ VLSs)]}2 ' (431) 



where Px is the power spectrum of the corresponding perturbation, X = 1Z for curvature perturbations and X = T 
for tensor perturbations, ji(x) are spherical Bessel functions 8] and 770 — 77^55 is the comoving distance between today 
and the last scattering surface (LSS) given by 

1 f 1 da , 

r)o - Vlss = 77- / /0 , o 1 i ' ( 4 - 32 ) 

-Ho 71 — vi2 r + il M a + S2a a 4 

where f2 r , J7m and Qa stand for the fraction of radiation, matter and cosmological constant in today's Universe. We 
find using zlss — 1100, 

3 - 296 (A ^ 

Vo - Vlss = -ft— • (4.33) 

-no 

Notice that k/Ho ~ dn / \hys(to) is the ratio between today's Hubble radius and the physical wavelength. The power 
spectrum for curvature (1Z) perturbations Pn(k) is given by eqs. (|4.21l) - (|4.24[) . 
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FIG. 9: The change AC e /C e on the CMB multipoles for £ = 1, . . . , 5. Upper plot: ACe/Ce vs. r - r* for < r - r* < 
0.2487963 . . .. Lower plot: ACe/Ce vs. r - r» for 0.0193 < To - t, < 0.2487963 .... r is the time when the BDic eq. (|477)l 
are imposed to the fluctuations. We choose to inside the fast-roll stage. ACe/Ce is positive for small to — r, and decreases 
with to becoming then negative. The CMB quadrupole observations indicate a large suppression thus indicating that 
to — t* > 0.05 ~ 10100 Tpianck- Our predictions here for the quadrupole and octupole suppressions are to be confronted 
with forthcoming CMB observations. It will be extremely interesting to measure the primordial dipole and compare with our 
predicted value. 



Inserting eq. (14.21)) into eq. (I4.31|) yields the Ci as the sum of two terms 



_ iBD AC, J,', x ' />;/■• .r; /;(x) <lr 



C^Cr + AQ , - 7 ^= J " " f : y , x = k{ m - VLSS ) = k/ K , (4.34) 
Ci J fi(x) dx 



where from eq.gSS]), k = H /3.296 . 



Mx) = x n °- 2 [j^x)] 1 . (4.35) 



and ji(x) stand for the spherical Bessel functions. 

The Cj s correspond to the standard BD power spectrum P^ D {k) eq. (|4.24p and the AC/ exhibit the effect of the 
transfer function D(k) on the C/. 

Using the transfer function D(k) obtained above eq. (14.221) . we computed the change on the CMB multipoles ACt/Ci 
for £ = 1, . . . , 5 as functions of the starting instant of the fluctuations tq. We plot ACg/Ci for 1 < £ < 5 vs. tq — t* 
in fig. [9l We see that ACi/Cg is positive for small tq — t* and decreases with To becoming then negative. The 
CMB quadrupole observations indicate a large suppression thus indicating that To — t* > 0.05 ~ 10100 Tpi anc k- 

Being D(k) < for low k as depicted in figs. [51 the primordial power at large scales is then suppresed and the low 
Ce decrease as seen from ea. (|4.34p . 
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ACe/Ci mainly originates from the peak of D(k) displayed in figs. [6] whose position moves to smaller k for 
decreasing t . Therefore, the primordial power suppression is less important for decreasing t and the CMB multipole 
suppression ACe/Ce less important as depicted in hgs. |U 

For small tq — t* < 0.05 the peak of D(k) grows significantly and ACg/Ci become positive, namely the low CMB 
multipoles are enhanced. 

It should be recalled that the observation of a low CMB quadrupole sparked many different proposals to explanain 
that suppression [l8j |. 

Besides finding a CMB quadrupole suppression in agreement with observations we provide here predictions 

for the dipole and octupole suppressions. Forthcoming CMB observations can provide better data to confront our 
quadrupole and octupole suppression predictions. It will be extremely interesting to measure the primordial dipole 
and compare with our predicted value. 

V. ANALYTIC FORMULAS FOR THE TRANSFER FUNCTION D(k). 

It is very important to dispose of analytic formulas for the transfer function D{k) in order to better understand the 
physical origin of its oscillations and properties as well as in the perspective of the MCMC data analysis. 

However, the mode equations (j2.24j) are not solvable in closed form for k ^ , not even for the approximated 
inflation solution eq. (|2.37[) which leads to the potential Vr(t) eq. (|2.43|) . 

The function D(k) must obey the general properties eq. (14.23[) . 

A. The primordial power spectrum vanishes for k — > and becomes the BD power spectrum for k — > oo 

The fluctuations equation (|2 . 24[) can be solved explicitly for k = 

s(v) = ci z(rj) + c 2 z(rj) J , (5.1) 

where c\ and c-i are arbitrary constants. 

The BDic eq. (|4.7l) introduce for k — >■ a l/\2k singularity in the mode functions. Thus, the mode functions must 
have the behaviour 

S K (k;r 1 ) k ^°^=[l + 0{k)] (5.2) 
V 2 fc 



where s{rf) is given by eq. (|5.II) . 

Inserting eq. (|5.2p into the BDic ea. (|4.7p yields for k — > 0, 

< \ 1 ds{rio) n 

s (Vo) = 1 , —, = , 

ar/ 

which determines the coefficients c\ and c 2 in eq. (|5.ip . We finally obtain 



and using eq. (|4.20p valid for rj — > 0~ when slow-roll applies 

hm S M = 1 . (5.4) 

>)^o- z(rj) z(r] ) 

The primordial power spectrum for k — > follows by inserting eq. (|5.2[) and eq. (|5.4|l into the general expression 
eq.flHES]), 



Pnik) ^ * lim 

\M PL J 2 it* n-to- 



z{rj) 



fe^o / m \ ( k 



M PL J \2n z(r) ) 
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We thus find in general that the power spectrum vanishes as fc 2 for k — > and therefore 

l + D(k) fc = O(F s+1 ) 

as stated in eq. (|4.23[) . This property is generally true except for the extreme slow-roll inflaton solution (sec. MI A[) 
with BDic imposed at 7/0 — — oo in which case D(k) vanishes identically for all k. 

For growing k the modes exit the horizon later on, during the slow-roll regime where eq. (|4.14ll applies. For large k 
the mode functions Sn as well as g„ TC behave as plane waves [eqs. (|4.6|) and (|4.18j) ] and therefore 

A n (k) = l , B n {k)=0 . Hence D{k) k . 



B. The transfer function D(k) when BDic are imposed during slow— roll. 

When the BDic eq. (|4.7j) are imposed during slow-roll at a finite time 770 we can use eq. (|4.14j) for the mode functions 
at r\ = t]q and we obtain, 

g— i k r/o 



t = A n {k) g„n(k;vo) + B n (k) 5* TC (fc;?7o) 

A K 



e -i k rjo 



-ik— =- = A n (k) g' v Jk; Vo ) + B n (k) g'^(k; Vo ) (5.5) 

V Z K 



which determines 



e — i k rjo g-i k r tQ 

A K(k)= , [g'*(k;r) )+ikg*(k;T] )] , B n (k) = - [g' v ^ (fc; t? ) + i k g„ n (fc; ry )] . (5.6) 

I v2 k i V2 k 

These coefficients satisfy eq. (|4.16|) and 

\A n {k)\ 2 + \B n (k)\ 2 = \ [\g' v Jk; Vo )\ 2 + k 2 \g^(k; Vo )\ 2 } 

Notice that the function g„(k; rf) eq. (|4.15p and the k factors in eq. (|5.6l) combine to produce functions A-jz(k) = An(k r/o) 
and Bn(k) = Bn(k t/q) that only depend on the product fc 770. 

We find from eqs. (|4.22p and (|5.6|) the corresponding transfer function which is a function of fc 770 too, 

l + D(k r /0 ) = 1 {|.g^(fc;ry )| 2 + fc 2 \g^(k;Vo)\ 2 ~ Re [z 3 " 2 ^ {g' 2 Jk; Vo ) + fc 2 9^(fc;»7o))] } (5.7) 

The functional dependence on fc 770 confirms the assertion in sec. HVBI that different initial times To lead to a rescaling 
in fc. 

In the fc 770 — > 00 limit two types of vanishing terms show up in D(k 770): (a) terms that strongly oscillate as 
e ±2% k no as they tend to zero and (b) non-oscillatory decreasing terms. Under integrals on fc, the terms of type (a) 
yield convergent expressions. We derive the non-oscillatory decreasing terms (b) by inserting the asymptotic behaviour 
of the Hankel functions eq. (|4.15|) Q in eq. (|5.7[) with the result 

D(k m) k =°° \ „ Hi + terms oscillating as e ±2t k 1)0 . (5.8) 
8 (fc 770) 4 

However, this approximation will not be valid for large enough fc since the modes at small enough wavelength will exit 
the horizon after the end of slow-roll where eq. (|5.8[) does not apply anymore. We recall that the occupation number 
\Bn(k)\ 2 (and therefore D(k)) must decrease faster than 1/fc 4 for fc — > 00 in order to ensure finite UV values for the 
expectation value of the energy-momentum fluctuations [2, 13]. 

The case Vfi — 3/2 is a good approximation which simplifies the expressions above. We obtain in this scale invariant 
case: 

• 1 -2ikrio 

An(k) = l + — --— 3 , B n (k) 



k 7 /0 2 fc 2 r)l ' V ; 2k 2 r]l 
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The transfer function is in this case, 

~. , cos 2x sin 2a; sin 2 a; „ ,„ , ,„ „, 

D(x) = — ^ 3- + ^^ . ^ = 3/2 , x = k m . (5.9) 

X X° X 

Eg. (|5.8p for v = 3/2 coincides with eq. (|5.9[) in the x — > oo limit, as it must be. 

Notice that the simple formula eq. (|5.9|) obeys the general properties eq. (|4.23j) . In particular, 

D(x) x =° -l + ^-x 2 + Q(x i ) . 
9 



VI. FIXING THE TOTAL NUMBER OF INFLATION E-FOLDS AND THE BOUND FROM ENTROPY 



It is very useful to plot the comoving scales of the cosmological fluctuation wavenumbers and the comoving Hubble 
radius together [see fig. HO] . One sees in this way how and when the cosmological fluctuations cross out and in 
the Hubble radius. The comoving Hubble radius is defined by Rh = l/[a(r) H(t)]. We display in Table VI the 
dependence of Rr on the scale factor a for all the relevant eras of the universe. 



Expansion stage 


Dependence of Rh on a 


Extreme Fast-roll 


a 2 


Fast-roll 


a 2 /%/a 6 + constant 


Slow-Roll inflation 


1/a 


Radiation Dominated 


a 


Matter Dominated 





TABLE V: Dependence of the comoving Hubble radius Rh = l/[a H] on the scale factor a for the relevant eras of the universe. 



The observed CMB quadrupole suppression can be easily explained if it exited the horizon by the end of fast-roll 
5, 6]. In that case, the modes which are horizon size today had wavenumbers &q ~ 11.5 m at horizon exit @. 
Combining this value of &q with the redshift since the pivot wavenumber exited the horizon, eqs. (|4.28l) . (|4.29l) and 
(|4.30|) . determines the total redshift since the beginning of inflation to be 



Ztot 



0.9 10 56 ~e 129 



Combining this value with the value of 1 + z r ~ 4 1 28 ~ e 66 by the end of inflation ea. (14.29[) yields a total number of 
Ntot — 63 inflation efolds. This value is very close to the minimal number of inflation efolds required to explain the 
entropy of the present universe due to photons and neutrinos 0): 



N tot > 62 .4 



Namely, this is the minimum number of inflation efolds compatible with the present entropy of the universe. 

In summary, assuming that the CMB quadrupole is suppressed because it exited the horizon by the end of fast-roll 
inflation fixes the total number of inflation efolds which turns to be 

N tot - 63 . 
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